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We show here a general model of phase separation in isotropic condensed matter, namely, a 
viscoelastic model. We propose that the bulk mechanical relaxation modulus that has so far been 
ignored in previous theories plays an important role in viscoelastic phase separation in addition to 
the shear relaxation modulus. In polymer solutions, for example, attractive interactions between 
polymers under a poor-solvent condition likely cause the transient gellike behavior, which makes 
both bulk and shear modes active. Although such attractive interactions between molecules of the 
same component exist universally in the two-phase region of a mixture, the stress arising from attrac- 
tive interactions is asymmetrically divided between the components only in dynamically asymmetric 
mixtures such as polymer solutions and colloidal suspensions. Thus, the interaction network between 
the slower components, which can store the elastic energy against its deformation through bulk and 
shear moduli, is formed. This unique feature originates from the difference in mobility between two 
components of a mixture. It is the bulk relaxation modulus associated with this interaction network 
that is primarily responsible for the appearance of the sponge structure peculiar to viscoelastic phase 
separation and the phase inversion: It suppresses the short-wavelength concentration fluctuations 
in the initial stage and causes the volume shrinking of a more viscoelastic phase. We also propose a 
simple general law of the stress division between the two components of a mixture, as a straightfor- 
ward extension of that obtained in polymer mixtures. We demonstrate that a viscoelastic model of 
phase separation including this new effect is a general model that can describe all types of isotropic 
phase separation including solid and fluid models as its special cases without any exception, if there 
is no coupling with additional order parameter. We show that this feature leads to a phenomenon 
of "order-parameter switching" during viscoelastic phase separation, even if it is driven by a single 
thermodynamic driving force. The physical origin of volume shrinking behavior during viscoelastic 
phase separation and the universality of the resulting spongelike structure are also discussed. 
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I. INTRODUCTION 



Phase-separation phenomena are widely observed in 
various kinds of condensed matter including metals, semi- 
conductors, simple liquids, and complex fluids such as 
polymers, surfactants, colloids, and biological materials. 
The phenomena play key roles in the pattern evolution 
of immiscible multi-component mixtures of any material. 
Thus, the phase-separation dynamics has been inten- 
sively studied in the past two decades from both the ex- 
perimental and the theoretical viewpoints From the 
concept of dynamic universality, phase-separation phe- 
nomena have been classified into various theoretical mod- 
els by Hohenberg and Halperin For example, phase 
separation in solids is known as "solid model (model B)" , 
while phase separation in fluids as "fluid model (model 
H)" 1^. For the former the local concentration can be 
changed only by material diffusion, while for the latter 
by both diffusion and ffow. It has been established that 
within each group the behavior is universal and does not 
depend on the details of material 

In all conventional theories of critical phenomena and 
phase separation, however, the same dynamics for the 
two components of a binary mixture, which we call "dy- 
namic symmetry" between the components, has been im- 



plicitly assumed [y^]- However, such an assumption of 
dynamic symmetry is hardly valid in various mixtures, es- 
pecially in a material group of "complex ffuids" . Recently 
we have found that in mixtures having intrinsic " dynamic 
asymmetry" between its components {e.g. polymer so- 
lutions composed of long chainlike molecules and simple 
liquid molecules), critical concentration fluctuation is not 
necessarily only the slow mode of the system and, thus, 
we have to consider the interplay between critical dy- 
namics and the slow dynamics of material itself . In 
addition to a solid and a fluid model, thus, we need the 
third general model for phase separation in condensed 
matter, which we call a "viscoelastic model" [^1). 

To describe the difference in elementary dynamics be- 
tween the two components of a mixture, we need a basic 
model that can treat the motion of each component sep- 
arately: such a model is known as a "two-fluid model". 
The basic dynamic equations of viscoelastic model have 
been derived to understand the coupling between the 
stress and diffusion and also the unusual shear ef- 
fects in polymer solutions ||9|-p^, which are known as 
"Reynolds effects", on the basis of a two-fluid model 
illjll- 

In this paper, we propose that we need some essen- 
tial modification to the "viscoelastic model" of phase 
separation described above: We believe [pjl that the 
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bulk relaxation modulus, which has been neglected (or, 
more strictly, not treated as an important physical fac- 
tor) in previous theories , plays an important role 
in viscoelastic phase separation as in gel phase separa- 
tion. This modified viscoelastic model can describe any 
kinds of phase separation in mixtures of isotropic con- 
densed matter without any exception, if there is no cou- 
pling with additional order parameter. In Sec. II, we 
review the theoretical derivation of a viscoelastic model 
based on a two-fluid model. In Sec. Ill, we discuss the 
internal modes of material itself and the resulting stress 
and how the stress is partitioned between the two com- 
ponents. In Sec. IV, we describe the basic equations 
of a viscoelastic model. In Sec. V, we discuss the ori- 
gins of asymmetric stress division, using a few examples. 
We also propose a general rule of stress division that is 
independent of material. In Sec. VI, we discuss the gen- 
erality of a viscoelastic model and demonstrate that all 
the models of phase separation in isotropic condensed 
matter are special cases of a viscoelastic model. In Sec. 
VII, we discuss the viscoelastic suppression of local con- 
centration fluctuations in polymer solutions, focusing on 
the roles of bulk relaxation modulus. In Sec. VIII, we 
demonstrate that characteristic features of viscoelastic 
phase separation can be explained by a simple concept of 
"order-parameter switching" , which originates from the 
general nature of the viscoelastic model. In Sec. IX, 
we discuss the universal nature of spongelike morphol- 
ogy characteristic of a dynamically asymmetric mixture 
containing a fluid as its component. In Sec. X, we con- 
clude our work. In Appendix, we briefly mention the 
applications of viscoelastic phase separation to material 
science. 



The volume average velocity v is given by 

(l}Vl + il~(l})v2. (2) 

The free energy of the system Fmix is given by 

Frm. = J dr[ficj,{r)) + ^{Vmfl (3) 

where /(</>) is the free energy per unit volume of a mix- 
ture with the concentration of the component 1. The 
form of /((/)) depends upon the system; for example, it is 
given by the Flory-Huggins free energy for polymer 
mixtures. Its time derivative can be written as 



9/ 
del) 



df 



-CV^<j>][V ■ {(l)vi)]dr 
(V • n) • vidr, (4) 



where V • 11 = '^V(|^ — CV'^cj)) and 11 is the osmotic 

tensor. We also assume here that the forces Fi acts on 
the component i. Thus, the Rayleighian to be minimized 
is 

R = j drl^p-^^^ + Icmvi - V2f + (V • n) • i/i 

-pV -v -vi- Fi-V2- F2]- (5) 

In the above, the term containing the pressure p is added 
to guarantee the incompressibility condition 



II. VISCOELASTIC MODEL BASED ON A 
TWO-FLUID MODEL 

We review here how a viscoelastic model can be de- 
rived on the basis of the two-fluid model [pT|-|l3| (see 
Ref. ||l^ on the details of the theoretical method). The 
model has originally be derived to describe the phase 
separation in polymers ||l2|,|l3|, but we believe that the 
viscoelastic model should describe phase separation in 
any dynamically asymmetric mixture, irrespective of the 
microscopic details of a system [^1). Thus, we focus here 
our special attention on how the most general version of 
the viscoelastic model can be derived. Let us consider 
a two-fluid model of a mixture of components 1 and 2. 
Let vi (r, t) and V2 [r, t) be the average velocities of com- 
ponents 1 and 2, respectively, and 4>{r, t) be the volume 
fraction of the component 1 at point r and time t. Here 
we assume the two components have the same density p 
for simplicity. Then, the conservation law gives 



V • z7 = 0. 



(6) 



dt 



-V-(M) = V-[(1 



(1) 



The condition that the functional derivatives of the 
Rayleighian with respect to vi and V2 be zero gives the 
following equations of motion: 

p^^^^ = -V ■Il-C{vi-V2) + <pVp + Fi, (7) 
p'-^^^^ = av.-V2) + {l^^)Vp + F2. (8) 
Thus, the average velocity v obeys 



dv - - -• 

Ptt = -V ■ n + Vp + Fi + Fa. 
at 



(9) 



In the quasi-stationary condition, the velocity difference 
between the two components, on the other hand, obeys 

vi-V2^ -i[(l - 0)V • n - (1 - <p)Fi + 0F2]. (10) 
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III. COUPLING OF INTERNAL MODES OF 
MATERIAL WITH DEFORMATION AND THE 
DIVISION OF THE RESULTING STRESS 

A. Origins of stress 

To obtain explicitly the form of the forces Fi, we need 
to understand how the stress is partitioned between the 
two components and also to have the microscopic expres- 
sion of the stress tensor of the material. The macroscopic 
total force F should be related to Fi and F2 as 



aivi + a2V2- 



(14) 



F 



F, +F2 



(11) 



Here <t is the total stress tensor, which is, in general, 
given by the constitutive equation of material. In a 
linear-response regime, the most general expression of cTy 
is formally written by introducing the time dependence 
of bulk and shear moduli in the theory of elasticity | |T9| ] 
as 



dt'[G{t-t')Kmt') 

) 



where 



dxi 



dv]. 



(12) 



(13) 



Here Vr is the velocity relevant to the rheological defor- 
mation and d is the spatial dimensionality. G{t) and 
K{t) are material functions, which we call the shear and 
bulk relaxation modulus, respectively. Here it should 
be noted that K{t) does not contain the bulk osmotic 
modulus, Kos — 0^(9^//(90^). We have the relation 
T] — Jp°° G{t)dt, where rj is the viscosity of the material. 

The second term of Eq. ^1^ ) is newly introduced to 
incorporate the effect of volume change into the stress 
tensor [^ , ^ . In a two-component mixture, the mode 
associated with V • tv can exist as far as Vr 7^ v, even if 
the system is incompressible. It should be stressed that 
its diagonal nature leads to the direct coupling with dif- 
fusion: note that the effective osmotic pressure is given 
by n-ff = (0§i - /) - dt'K{t ~t')W- Vrit')- We be- 
lieve, thus, that this term is important even in the case of 
polymer solution, as described later, although this term 
has so far been ignored (or, more strictly, not treated 
as an important physical factor) in the previous theories 



B. Estimation of the stress division parameter au 

Here we consider the physical meaning of tv. In a 
linear-response regime, the rheological velocity iv is gen- 
erally given by the linear combination of vi and V2 |l2,n : 



Then, the next problem is how the stress is partitioned 
between the two components. Since F-Vr should be equal 
to Fi ■VI+F2-V2 in the Rayleighian, we have the following 
stress division: 



Fi = aiF, F2 — 0L2F. 



(15) 



Here ai + a2 = 1 from Eq. (11) 



C. Direct Estimation of Fk 

Here we consider the meaning of the forces from a dif- 
ferent viewpoint. The forces acting on the component 
k are (i) the friction between the component k and the 
other component due to their relative motion and (ii) 
the rheological coupling between the component k and 
the surrounding rheological environment including the 
component k itself. This can be easily understood by 
considering gel that is composed of polymer network and 
solvent, as an example: The motion of polymer is affected 
by the two forces, namely, the friction force against sol- 
vent and the elastic force due to the network deformation. 
Thus, it is natural to think that Fk (/c = 1, 2) corresponds 
to the force of type (ii), namely, the force acting on the 
component k by the motion of the component k (vk) it- 
self and not by that of the other component. Thus, we 
assume that Fk is linear in Vk'. 



where 



dt'[GW(i-i')Kfe'(i') 



+K^'^\t-t')iV ■Vkit'))S,j], 



(16) 



(17) 



(18) 



Here the unknown factors become the functional shapes 
of G^'^^t) and K''^\t) for the motion of the component 
fc, instead of ak- 

Using the stress division parameters au , we obtain the 
following relations between G and G'-'^^ and also between 
K and K^^^: 



K = 



G(i)G(2) 



k2G(2) ' 



X(l)if(2) 



(19) 
(20) 



These functions G^^) and X^^) express the rheological 
properties of the material responding to the velocity Vk- 
The fact that vi and V2 are coupled with each other 
makes the estimation of the rheological functions diffi- 
cult. 
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D. Qfe or Fk 

The estimation of the stress division of F into Fi and 
F2, which is given by the stress division parameter ak, 
formaUy looks simpler than the direct estimation of the 
above rheological functions, G'-'^-' and X*^'^'; this is true 
for a nearly symmetric stress division, as in the case 
of polymer mixtures: More strictly speaking, when the 
dynamics of both components is governed by the same 
mechanism, ak is useful to estimate the stress division. 

However, the physical meaning of the latter is clearer 
than the former in a case of strongly asymmetric stress 
division: In polymer solutions, for example, the mech- 
anism of polymer dynamics is essentially different from 
that of solvent dynamics; and, thus, the rheological func- 
tions G'-'^-' and K^''^ are more useful than ak (see Sec. 
VI). 



IV. BASIC EQUATIONS OF A VISCOELASTIC 
MODEL 

Here we summarize the basic equations describing a 
viscoelastic model: 



dt 

Vi - V2 

dv 



-v-(0t/)-v-[<^(i-(/.)(t/i-t;2)], (21) 
1 



-[v-n- v-(T 



(1) 



p 



dt 



C ' ' 1- 

-V • n + Vp + V • cr(i) + V • (t(2 



V-cr(2)], (22) 
(23) 



We also need Eq. (|), Eq. Eq. (|T^), and the infor- 
mation on the stress division. Here it should be noted 
that we need the phenomenological or microscopic the- 
ories describing the forms of G{t) and K{t) or those of 
GC^) and K<'^\ Since we derive the above basic equa- 
tions, relying only upon a two-fluid model, they should 
be independent of types of material and quite general. 



V. ASYMMETRIC STRESS DIVISION 

Here we first focus on some specific problems to have 
more deep insight on the origins of stress and how the re- 
sulting stress is partitioned between the two components 
of a mixture. Then, we consider some general features of 
asymmetric stress division. 



A. Examples of asymmetric stress division 

1. Polymer solutions: cases of good and 6 solvents 

In this case, the stress division has also already been 
given in literature |ll-lj]: ai = 1 and 02 = 0, provided 



that the component 1 is a polymer and 2 is a solvent. 
However, this division itself is based indirectly on the es- 
timation of G^'^-' since the stress division parameters can- 
not be determined precisely. It should be stressed that 
there is no direct way to determine ai and a2 from the 
first principle. Thus, there is no firm basis for Vr — vi in 
polymer solution, as pointed out by Doi and Onuki Jl2[ , 
although it looks natural physically. In this case, we be- 
lieve that Eq. (|l^) is more useful and physically easier to 
understand than Eq. (p^), as described in Sec. HI. D. 

The topological entanglement can be felt only by poly- 
mers, and never by solvent molecules. Then G*^^-'(t) is ap- 
proximately given by the existing polymer-solution the- 
ory H [G(i)(0 ^ G{t)] and K^^\t) - IMI, if the 
solvent is not poor. On the other hand, G^'^'> (t) IS an 
extremely fast decay function and, thus, the solvent vis- 
cosity r]2 is obtained as 772 = dtG^^\t). Since the 
solvent has no large internal degrees of freedom, we can 
safely assume that K^'^'' = 0. 



2. Polymer solutions: a case of poor solvent 

It should be stressed that phase separation of poly- 
mer solution always occurs in a poor solvent; and, thus, 
the case of poor solvent is extremely important when we 
consider critical phenomena and phase-separation phe- 
nomena. Unfortunately, however, there do not exist any 
established theories that describe quantitatively the poly- 
mer dynamics in a poor solvent. In a poor solvent, we 
need to consider the attractive interactions between poly- 
mer chains seriously. Thus, the most natural model is a 
transient gel model |^,^ in which the interpolymer at- 
tractive interaction produces the temporal contact point 
(crosslinking) between polymer chains. If we assume that 
the life time of the temporal contact between chains is , 
we expect that the bulk relaxational modulus K^^'^ (t) has 
a relaxation time of the order of r^,. likely obeys the 
Arrhcnius-type law: = t!^ exp{E/kBT) {kB- the Bolz- 
mann constant). Here the bonding energy E is likely 
proportional to Tg — T (Tg: the 6 temperature) near Tg, 
and also dependent upon the distance and orientation of 
the relevant segments. Even in polymer solutions, thus, 
we expect that K^^^{t) plays an important role in vis- 
coelastic phase separation in contrast to the previous the- 
ories |]lT|-p^ : The transient network of topological origin 
itself ( entanglement effects ) might not lead to the hulk 
relaxation mode, while the transient network formed by 
attractive interactions does clearly makes the bulk relax- 
ation mode active. We believe that it is this mode that 
is primarily responsible for the volume shrinking behav- 
ior of a more viscoelastic phase observed in our experi- 
ments [|[^]. The characteristic decay time of G'^'(<) is 
likely longer than Tx since the reptation-like motion is 
additionally required for the polymer motion under the 
temporal network formed by the attractive interaction. 
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We need further theoretical studies to have quantitative 
expressions of G{t) and K{t) in a poor solvent. 

Relating to the above transient gel model, we speculate 
that a polymer solution behaves as physical gel univer- 
sally at least at a high polymer concentration under a 
strongly poor-solvent condition. The existence of special 
junction points is not a prerequisite to the formation of 
such a transient gel. The transient pairing of any parts 
of two chains can be regarded by a temporal crosslink- 
ing. Any pair of segments of polymer chains can form 
a temporal crosslinking point, irrespective of inter or in- 
trachain. The probability of its formation is determined 
by the balance between the intersegment attractive in- 
teraction depending upon the geometrical configuration 
of chains and the thermal energy. Thus, the most prob- 
able candidate of the contact point is an entanglement 
point. The sol-gel transition can be given simply by the 
criterion that the transient network formed by interpoly- 
mer attractive interaction is percolated at any moment. 
A polymer chain having at least two contact points with 
other different polymers plays a role as a junction point of 
usual physical gel. If we assume simply that the topologi- 
cal entanglement point is the only candidate for a tempo- 
ral crosslinking point, the criterion for physical gelation 
is given by 



^exp(£;/fcsT) >2, 



(24) 



where N is the degree of polymerization of polymer and 
Nf, is the degree of polymerization between entanglement 
points. Further quantitative studies along the above line 
are highly desirable. 



3. Polymer mixtures 

In the case of a mixture of polymers 1 and 2, whose de- 
grees of polymerization are Ni and respectively, the 
stress produced by the motion of polymer 1 can be differ- 
ent from that by the motion of polymer 2. Intuitively, the 
motion of a longer chain causes stronger stress than that 
of a shorter chain does. Following the Brochard theory 
on mutual polymer diffusion [ p3| , which is based on the 
reptation theory that mainly deals with the effect of topo- 
logical constraints (tube) on entangled polymer chains, 
Doi and Onuki have explained how the stress should 
be divided by the two polymers with different lengths. 
According to them. 



Vr 



Oik 



vt — aivi + a2V2, 

Cfc ^ 0iVfc 
C1+C2 ^Ni + {1-^)N2 



(25) 
(26) 



fraction of the component fc, Co is the microscopic friction 
constant, and is the average degree of polymerization 
between the entanglement points. The resulting stress 
division is given by Ffc = a^V ■ cr Here it should be 
stressed that Cfc is the very essential quantity in the sense 
that it represents the coupling strength between the com- 
ponent k ( the volume fraction of 4>k ) and the surrounding 
rheological environment. 

Near and below a critical point Tc, however, we also 
have to consider the role of the attractive interactions 
between the same kind of polymers, which increases the 
rheological coupling, namely, Cfc- For example, this leads 
to the slower diffusion constant than that predicted by 
a reptation theory which concerns only topological ef- 
fects and neglects energetic interactions between poly- 
mers. The inclusion of energetic interactions is a prereq- 
uisite to the more precise description of polymer dynam- 
ics. However, the following fact should be mentioned: 
High-molecular weight polymer mixtures often mix at a 
lower temperature, and demix at a higher temperature. 
In such a case, the energetic interactions likely play more 
important roles in polymer dynamics in the one-phase 
region rather than in the dcmixing region. 



4- A mixture of components having very different 
glass-transition temperatures (Tg) 



In this case, we also expect an asymmetric stress di- 
vision since the two kinds of the component molecules 
are expected to feel very differently the rheological envi- 
ronment as in the case of polymer mixtures, even if the 
mean-field rheological environment surrounding them is 
the same. We have actually observed viscoelastic phase 
separation in a mixture of polymers having very differ- 
ent Tg, whose behavior is essentially the same as that of 
polymer solutions |^ . It is easy to imagine that a high Tg 
component has less friction with the local rheological en- 
vironment than a low Tg component. Recent theoretical 
studies on supercooled binary liquids HJ] based on the 
mode-coupling approximation support such a picture. If 
we introduce formally the coupling strength Cfc for the 
component k that is proportional to (f>k, the stress di- 
vision can be expressed by the same relation as in the 
above case of polymer mixtures. The mean-field rheolog- 
ical environment in this problem of glass transition is the 
The concept of "cage" in 



where vt is the tube velocity. Here Cfc is the friction 
of the component k with the tube surrounding it. Cfc is 
given by Cfc = (j)k{NkC,[)/Ne) where (j)k is the volume 



so-called "cage" |25 . The concept of "cage" in glass tran- 
sition is quite similar to the concept of "tube" in polymer 
mixtures. The escape time of a molecule from "cage" or 
"tube" gives the relaxation time of G{t) in both cases. 
Unfortunately, we do not have the reliable theoretical 
basis even for a simple liquid-glass transition; and, thus, 
it is difficult to have specific quantitative expressions for 
Cfc at present. Phenomenologically, however, it is known 
that G{t) = Goexp[-(i/r)^] andr ~ tq exp(B/(r-ro)), 
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where (3 is the stretching parameter (0 < /3 < 1) and Tq 
is the so-called Vogel-Fulcher temperature. 

It should be stressed again that we have to take into 
account the effects of attractive interactions between the 
same species on their dynamics below Tc pq| . 



5. Colloidal suspensions 

It is well known that the addition of enough non- 
absorbing polymer to an otherwise stable colloidal sus- 
pension can induce phase separation via the depletion 
mechanism. Colloidal suspensions form a transient gel 
state in the initial stage of phase separation p7|] in much 
the same way as polymer solutions do. We believe that 
the essential features of colloid phase separation can also 
be well described by our viscoelastic model. The dynamic 
asymmetry in colloidal suspensions simply comes from 
the size difference between colloids and solvent molecules. 



B. Physical origin of asymmetric stress division 

Here we consider a problem of what is the most basic 
physical factor that is responsible for asymmetric stress 
division, on the basis of intermolecular or interparticle 
interactions. The network of attractive interaction is 
universally formed when a mixture is quenched into its 
metastable or unstable state since there exist attractive 
interactions between the same components. In dynam- 
ically symmetric mixtures, the interaction network al- 
ways relaxes in its equilibrium state much faster than 
the phase-separation process. In dynamically asymmet- 
ric mixtures, however, the relaxation time of the inter- 
action network is different between the two components 
because of the mobility difference. This consideration 
based on microscopic interactions leads to the conclusion 
that the dynamic asymmetry between the components of 
a mixture is the essential origin of asymmetric stress di- 
vision. Thus, the phase-separation behavior of any dy- 
namically asymmetric mixtures including the above 1-5 
should be essentially the same and described by Eq. 



given by {vr—^k) , where Cfc is the average friction of the 
component k and the mean-field rheological environment 
at point r, where the volume fraction of k component 
is (j)k{r). Here Cfc = (/"feC™ and is proportional to the 
friction between an individual molecule of the component 
k and the mean-field rheological environment, which we 
call the generalized friction parameter. Because of the 
physical definition of the mean-field rheological environ- 
ment, the two friction forces should be balanced. This 
fact guarantees that the rheological properties can be de- 
scribed only by Vr as in Eq. (|l|). Thus, we have the 
following relation, in general: 



Cl(Wr - Vi) + (2{Vr - V2) = 0. 



(27) 



From Eqs. (|l|) and we obtain the general expres- 
sion of the stress division parameter ak '■ 



(28) 



The above relation is consistent with a simple physical 
picture that the friction is only the origin of the coupling 
between the motion of the component molecules and the 
rheological medium. The above relation is a straightfor- 
ward extension of the stress division in polymer mixtures 
Jl2t , where Vr is the tube velocity vt- We expect that 
this relation holds, irrespective of the microscopic de- 
tails of rheological models, and, thus, we can apply it to 
a mixture of any material, the motion of both of whose 
components is described by a common mechanism. How- 
ever, it should be stressed that this relation is not useful 
for mixtures whose components have essentially different 
mechanisms of molecular motion as in the case of poly- 
mer solutions: In polymer solutions, for example, the 
motion of polymers is essentially different from that of 
solvent molecules in the mechanism. In the standard 
rheological theory of polymer solutions, G{t) itself is es- 
timated as the sum of the polymer contribution and the 
solvent contribution. Thus, the stress division cannot be 
simply described by Eq. (psj). In such cases, Eq. ( |l7| ) is 
more useful, as mentioned in Sec. HI. D. 



C. A general rule of stress division 

On the basis of the above examples, we discuss a gen- 
eral rule of the stress division in viscoelastic matter. Here 
we do not consider elastic matter where the elastic cou- 
pling also plays an important role in addition to the fric- 
tion. In the preceding discussion, we obtain the gen- 
eral relation given by [Eq. (p^)], Vr = aivi -\- a2V2, with 
ai + Q!2 = 1. For the relative motion of the component 
k having the velocity of Vk to the mean-field rheological 
environment having the velocity of tv, the friction force is 



VI. GENERALITY OF A VISCOELASTIC MODEL 



Next we briefly discuss the generality of the above vis- 
coelastic model described by Eqs. (pT|), (|2^), and ( [2^ ) 
pO{ . This model including the bulk volume relaxation 
mode is quite a general model, as shown below. We de- 
scribe below how the viscoelastic model reduces to vari- 
ous models under some assumptions. 
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A. Elastic solid model 

If we assume that G{t) — ^{(p) {j-r. shear modulus) and 
K{t) — Ki,{ip) {Kb: bulk modulus) and t/ = 0, this model 
reduces to the model of elastic solid model ||2^ . Since the 
time integration of the velocity becomes the deformation 
u, the stress is given by 



duj 
dxj 



'-[V ■u)5,, 



(29) 



Thus, the basic kinetic equation is given by 



dt 



= V • 



1-^ 



-[v-n- v-cr(i) 

+ -^V-rT^% (30) 



In this case, the softer phase form a networklike phase be- 
cause the deformation of the softer phase costs less energy 
than that of the harder phase |2|]. It should be stressed 
that the force balance condition plays no roles in deter- 
mining the morphology. This fact causes the striking dif- 
ference in morphology between an elastic solid model and 
an elastic gel or asymmetric viscoelastic model Q : In the 
former the softer phase forms the networklike structure, 
while in the latter the harder phase does. 



B. Solid model 

If we assume the dynamic symmetry (no dependence 
of /i and Ki, on 0) further, it reduces the solid model 
(model B j2]). This is because we have the symmetric 
stress division as (1 — (t))Fi — <j)F2. Here it should be 
noted that the condition ^ = Ki, = is unnecessary and 
only the symmetry in elastic properties between the two 
components is required. In this case, the basic equation 
becomes the simplest diffusion equation: 



1- 



-[v-n] 



C. Symmetric Viscoelastic model 



(31) 



bulk relaxation modulus under the incompressibility con- 
dition (V • t/ = 0). The basic kinetic equations are given 

by 



(1 



dv ^ 

p— ^ -v-n + Vp + v-cr. 

at 



-V n. 



(32) 
(33) 



Since Vr = v, the gross variables describing the dynamics 
are only cj) and v. Here it should be stressed that the rhe- 
ological function G{t) does not depend upon the location 
r because of the dynamic symmetry. Using the relation 
V • u = 0, thus, 



dt'G{t - t')V^v{t') 



We have also the relation 



Vl 



V2 = 



C 



-[v-n]. 



(34) 



(35) 



although it is unnecessary for solving the problem. 

This model describes the dynamics of dynamically 
symmetric polymer mixtures. It should be stressed that 
this model is different from a fluid model (model H) de- 
scribed below. Thus, there remains a possibility that 
there is a new polymer effect relating to this model: For 
example, we expect an unusual feature in the initial stage 
of phase separation where the deformation rate is large. 
The polymer effect corresponding to this was first pointed 
out by de Gennes ||2^ , and has very recently been studied 
in detail by Kumaran and Frederickson pc|]. 



D. Fluid Model 

If we assume that the deformation is much slower than 
the internal rheological time of the material for the above 
model, we further have the relation V • cr = rfS/^v. Thus, 
the model reduces to the fluid model (model H |^): 

d(h , , ^ (hil- <bf ^ 

Ji^-V.{^v) + \7- V ■ n, (36) 



dv ^ o 
p— ^ -V ■U + Vp + ri^^v. 
at 



(37) 



If we assume only the dynamic symmetry between the 
two components of a mixture, it reduces to a new "sym- 
metric viscoelastic model" . In this case, we have a trivial 
stress division: Fi = (f>V ■ cr and F2 = (1 — 0)V • cr. 
Namely, ai = and a2 — 1 — (j). The rheological 
functions G^'^-' can be estimated as G'-^-'(i) = 4'G{t) and 
G'^'^\t) = (1 — (j))G{t). In this particular case, Vr = 
and, accordingly, there should be no contribution of the 



E. Elastic gel model 

If we assume only G = p,{(j)) and K — Ki,{(f)), it reduces 
to the elastic gel model [31 3^ that describes phase sep- 
aration in elastic gel. The basic equations are essentially 
the same as those of viscoelastic phase separation in poly- 
mer solution [Eqs. (pl|), (p^), and (|2^)], except that the 
stress tensor is given by 
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.(2) 



'72 I 



,dv2 



In this case, when the elastic energy overcomes the mix- 
ing free energy, phase separation stops and the coarsening 
of domains is pinned. 



F. Generality and intrinsic nonuniversality 



Since any phase separation in aU isotropic condensed 
matter can be classified into solid, elastic solid, elastic gel, 
symmetric and asymmetric viscoelastic, and fluid mod- 
els, the above viscoelastic model (see Sec. IV) includ- 
ing both shear and bulk relaxation stresses should be a 
universal model describing phase separation and critical 
phenomena in isotropic matter without any exception. 

However, this model is not universal in the usual sense 
of critical phenomena since it requires some microscopic 
theories describing the rheological properties of the mat- 
ter. In the extreme limit of the strong dynamic asym- 
metry, the elementary slow dynamics (internal mode) of 
material affects the critical fluctuation even near the crit- 
ical point, in contrast to the concept of the dynamic uni- 
versality 1^. Thus, there is a possibility that we cannot 
experimentally approach to a critical regime where the 
order parameter dynamics is only the slow mode in the 
system @|. 

The viscoelastic effects can be parameterized by the 
so-called viscoelastic length ^ [D^TtY^"^ originally 
introduced by Brochard and de Gennes H], where Tt is 
the characteristic time of rheological relaxation and 
is the diffusion constant. For the length scale longer 
than ^i,e, concentration fluctuations decay by diffusion, 
while in the length scale shorter than viscoelastic effects 
dominate [0jl2|jl|]. This is a simple mapping of the dy- 
namic crossover that in the time scale longer than Tt dif- 
fusion dominates concentration fluctuations while in the 
time scale shorter than Tt viscoelastic effects dominate. 
The critical regime is, thus, described by the condition 
f ^ £,ve, where ^ is the correlation length of concentra- 
tion fluctuations. This condition can also be written as 

^ Tt, where is the characteristic time of the criti- 
cal concentration fluctuations. Thus, we need to consider 
whether we can easily approach to the critical regime that 
is defined by the above criterion, in a viscoelastic system. 
The more detailed consideration on this problem of in- 
trinsic nonuniversality will be described elsewhere. 



VII. VISCOELASTIC SUPPRESSION OF LOCAL 
CONCENTRATION FLUCTUATIONS IN 
POLYMER SOLUTIONS 

A. Fluctuation suppression due to bulk relaxation 
modulus 

According to the continuity equation, Eq. (I), we have 
the following relation. 



— = -V0 • t7i - (/>V • vi. 
ot 



(38) 



Here we assume that the component 1 is a polymer and 
the component 2 is a solvent. In the above equation, 
the first term of the right-hand side simply describes the 
translational transport of polymers to a point f by the 
locally uniform velocity field vi , while the second one de- 
scribes the polymer diffusion toward or from a point r. 
Thus, the first term is associated only with the change 
in the spatial pattern of the concentration distribution, 
while the second term is responsible for the change in 
the concentration distribution itself. In the initial stage 
of phase separation, the major process is the diffusion 
process leading to the change in the concentration dis- 
tribution itself and there are few changes in the spatial 
pattern of the concentration distribution. Neglecting the 
first term in Eq. (B8), thus, we have the relation 



dt 



/(/) ^ -V • VI. 



(39) 



The left-hand side of the above equation is inversely pro- 
portional to the characteristic time of the concentration 
change, r^. On the other hand, the bulk relaxation mod- 
ulus K{t) that is directly coupled with V • vi has a char- 
acteristic decay time of t^ , which is related to the charac- 
teristic time of the transient crosslinking between poly- 
mer chains in a poor-solvent condition. If t^ ^ t^, the 
rapid growth of concentration fluctuations characteristic 
of spinodal decomposition is suppressed severely and may 
even be prohibited. If t^ <C t^, on the other hand, there 
are few elastic effects and the concentration fluctuations 
can grow as in usual spinodal decomposition. 



B. Spinodal decomposition vs. nucleation and growth 

If the above mechanism of the suppression of concen- 
tration fluctuations works efficiently, spinodal decompo- 
sition cannot proceed further after the crossover between 
Tj; and r^. Thus, the type of phase separation switches 
from spinodal decomposition to nucleation and growth. 
More strictly speaking, the concentration fluctuations 
grow only locally in the solvent-rich region where t^ is 
short and the diffusion is easy to take place. Such be- 
havior is actually observed in our experiments The 
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crossover from an initial fluid state to a transient gel state 
likely takes place almost immediately after the tempera- 
ture quench, since the long-range diffusion or motion of 
molecules are not required for the formation of a tran- 
sient network (see also Sec. VIII. C). After the for- 
mation of a transient gel state, the "mechanical insta- 
bility" that is a universal feature of the "soft" network 
of attractive interactions leads to the nucleation-growth- 
like behavior. In the following, we discuss the concept of 
"order-parameter switching" resulting from the crossover 
between the characteristic phase-separation time and the 
internal rheological time. 



For Td ^ Tfs the primary order parameter is the com- 
position in usual classical fluids, while for < Tts it is 
the deformation tensor [dij = duj/dxi -I- dui/dxj) as in 
elastic gels. In the elastic regime, the force terms can be 
included in the Hamiltonian as in the case of gel. Then 
the free energy functional is formally written only by the 
deformation tensor dij as f{dij). Thus, we can say that 
the order-parameter switching is a result of the compe- 
tition between two time scales characterizing domain de- 
formation Td and the rheological properties of domains 
Tts- This is a kind of viscoelastic relaxation in pattern 
evolution. 



VIII. ORDER-PARAMETER SWITCHING 

A. Order-parameter switching between composition 
and deformation tensor 

Here we consider the dynamic process of viscoelastic 
phase separation on the basis of the viscoelastic relax- 
ation phenomena described by Eq. (]T^). The quanti- 
tative feature of the dynamics can be understood on 
the basis of a concept of "order-parameter switching" 
pof . Phase separation is usually driven by a thermo- 
dynamic force and the resulting ordering process can be 
described by the temporal evolution of the relevant or- 
der parameter associated with the thermodynamic driv- 
ing force. The primary order parameter describing phase 
separation of a binary mixture is a composition differ- 
ence between the two phases. Besides exceptional cases 
where phase separation and other ordering processes 
such as superfluidization, gelation, liquid-crystallization, 
and crystallization, simultaneously proceed [in other 
words, there are more than two kinds of thermodynamic 
forces (order parameters)], a phase-separation process is 
usually characterized by a single order parameter. In 
the viscoelastic model, on the other hand, the phase- 
separation mode can be switched between "fluid mode" 
and "elastic gel mode". This switching is likely caused 
by the change in the coupling between stress fields and 
velocity fields, which is described by Eq. (|l^): Equa- 
tion (^2|) tells us that these two ultimate cases, namely, 
(i) fiuid model (k^^-, V • w,- ~ const) and (ii) elastic gel 
model {G{t), K{t) ^ const), correspond to Tts ^ and 
Tts *C Td, respectively. Here Td is the characteristic time 
of deformation, and Tts is the characteristic rheological 
time of the slower phase. There can be two types of Tts ■ 
One is associated with the characteristic decay time of 
G{t), while the other with that of K{t). We think that 
the former is generally longer than the latter in polymer 
solutions, as mentioned in Sec. V. A. 2. As described in 
the preceding section, we believe that the latter plays an 
important role in the suppression of concentration fluc- 
tuations in the initial stage. 



B. How does the order-parameter switching occur? 

We next consider how Tts and Td change with time 
during phase separation. In the initial stage, the ve- 
locity fields grow as v ^ {kBTK/3ri£^)A(fP , where A(j) 
is the composition difference between the two phases, 
and ^ is the correlation length, or the interface thick- 
ness. Since A(f) approaches to 2(f)e [(pf. :the equilibrium 
composition) with time, this expression of v reduces to 
the well-known relation v ^ ^j-x] (7: interface tension) 
in the late stage [note that 7 - k^TK(2(\)e)^ ji^. Thus, 
the characteristic deformation time Td changes with time 
as Td ^ R{t)/v{t) - R{t) / Acj){t)'^ . In the initial stage, 
the domain size does not grow so much with time while 
A(j3 rapidly increases with time; and, accordingly, Td de- 
creases rapidly. On the other hand, Tts increases steeply 
with an increase in A0, reflecting the increase in the 
polymer concentration in a polymer-rich domain. Thus, 
Tts becomes comparable to Td in the intermediate stage of 
phase separation. Once Td < Tts , the slower phase cannot 
follow the deformation speed and behaves as an elastic 
body: The elastic energy dominates the coarsening pro- 
cess in the intermediate stage. Next we consider the late 
stage. Since A(j) approaches to 20e and becomes almost 
constant in the late stage, Td (~ Rr]/j) increases with an 
increase in R while Tts becomes almost constant. Thus, 
Td becomes longer than Tts again. In short, Td ^ Tts 
in the initial stage, Td < Tts in the intermediate stage, 
and Td 3> Tts in the late stage again. Accordingly, the 
order parameter switches from the composition to the 
deformation tensor, and then switches back to the com- 
position again. This is the first case of "order-parameter 
switching" during an ordering process driven by a single 
thermodynamic driving force, to our knowledge. 

C. Further consideration on G{t) and K{t): Kinetics 

of transient gel formation 

Here we briefiy consider viscoelastic functions, G{t) 
and K{t), which are dependent upon material. In the 
above discussion, Tts is assumed to be a function of A^. 
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However, this picture is not necessary true. In poly- 
mer solutions and colloidal suspensions, for example, the 
transient gellike structure is likely formed very quickly 
after the quench. This is because the diffusion of poly- 
mers or colloidal particles over a large length scale is not 
required to form the interaction network. The diffusion 
length scale / is the order of polymer or particle size, a, 
near the critical composition {(j)c)', and, thus, the time 
required to form network is a^/Da, where Da is the 
diffusion constant of a polymer or particle. In such a case, 
Tts very rapidly increases to the order of within a time 
of ^ a? /Da after the quench. In such a case, the first 
order parameter switching from composition to deforma- 
tion tensor occurs within a very short period (~ a? /Da) 
after the quench: The system enters into an elastic regime 
just after the quench. 

The diffusion length scale I increases with a decrease in 
polymer or colloid concentration 0. If a percolated net- 
work cannot be formed within a sufficiently short time, 
a quasi-homogeneous transient gel state is not realized 
due to phase separation, and, thus, a networklike phase- 
separated pattern is not formed; instead, a droplet pat- 
tern is formed |^]. This criterion may give the threshold 
composition between droplet and network phase separa- 
tion. 

Finally, it should be stressed that to describe this net- 
work formation process and the resulting change in vis- 
coelastic functions, G{t) and K{t)^ we need a microscopic 
model of each system. 



Lhoie- Neglecting the time dependence of Lhoie, we get 
a very rough estimation of the volume shrinking time 
as Tshrink ^ ^hoie/^^ which Can be a very long time. 
The estimation of Lhoie including its time dependence re- 
quires the stability analysis under the influence of phase 
separation. We need further study on this problem. 



E. Pattern selection: elastic energy vs. interface 
energy 

Since the deformation tensor intrinsically has geomet- 
rical nature, the pattern in the elastic regime is essen- 
tially different from that of usual phase separation in 
fluid mixtures. The domain shape during viscoelastic 
phase separation is determined by which of elastic and 
interface energy is more dominant. Roughly, the elastic 
energy is estimated as fie^R'^ (e: strain and d: spatial 
dimensions) for a domain of size R, since it is the bulk 
energy. On the other hand, the interface energy is es- 
timated as jR'^^^. For macroscopic domains, thus, the 
elastic energy is much more important than the interface 
energy in the intermediate stage where < Tts- Ac- 
cordingly, the domain shape is determined by the elastic 
force balance condition (V ■ cr^"' ~ 0), which leads to 
networklike or spongelike morphology. In the initial and 
late stages of phase separation where ^ Tts, on the 
other hand, the interface energy dominates the domain 
shape since ^ ^ 0. 



D. Volume shrinking behavior: Absence of 
self-similar pattern growth in viscoelastic phase 
separation 

Because of the order-parameter switching, there is no 
self-similarity in pattern evolution of viscoelastic phase 
separation. In the elastic regime, further, the volume ra- 
tio between the two phases changes with time 
and, thus, there is no proportionality between interdo- 
main distance and domain size. This behavior even leads 
to the phase inversion when the more viscoelastic phase 
is a slightly minority phase in equilibrium: in the initial 
stage a less viscoelastic phase forms droplets, while in the 
final stage a more viscoelastic phase does. This means 
that there are at least two length scales that change dif- 
ferently with time. This is also related to the absence of 
any characteristic length scale in elastic deformation. 

Here we make a rough estimate of volume shrinking ki- 
netics. The elastic regime should be analogous to the vol- 
ume shrinking of gel. The characteristic shrinking time 
Tshrink IS likely dependent upon the characteristic length 
scale L as / D {D: gel diffusion constant |^). It should 
be noted that only the length scale in this problem is the 
length scale associated with the "mechanical instability" , 
namely, the characteristic distance between solvent holes, 



IX. UNIVERSALITY OF SPONGE 
MORPHOLOGY CHARACTERISTIC OF A 
MIXTURE HAVING BULK COMPRESSION 
MODES 

A. Universality of sponge morphology 

We discuss here the universal nature of a spongelike 
morphology (or the formation of a continuous structure 
by a minority phase) and its physical origin. It is known 
that gel undergoing volume-shrinking phase transition 
forms a bubble- like structure [^ , ^3|j3^ ] . The competition 
between phase separation and gelation or glass transition 
also causes spongelike morphology p5|-^. The physi- 
cal origin of the appearance of a honeycomb structure 
in plastic foams {e.g. polystyrene foam and urethane 
foam) is also similar to ours. All these processes have 
a few common features: (i) A mixture contains a fluid 
as a component, (ii) Holes of a less viscoelastic fluid 
phase (gas in plastic foam, water in gel, solvent in poly- 
mer solution, and so on) are nucleated to minimize the 
elastic energy associated with the formation of a hetero- 
geneous structure in an elastic medium, (iii) Then, a 
more viscoelastic phase decreases its volume with time. 
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This volume shrinking process is dominated by the trans- 
fer (diffusion or flow) of a more mobile component under 
stress fields, from a more viscoelastic phase to a less vis- 
coelastic phase. The above picture suggests the possibil- 
ity that a spongelike structure is the universal morphol- 
ogy for phase separation in systems in which only one 
component asymmetrically has elasticity stemming from 
either topological connectivity or attractive interaction. 

We also point out |^ the similarity of these patterns in 
condensed matter to the spongelike structure of the uni- 
verse (the large-scale galaxy distribution) . We spec- 
ulate that the gravitational attractive interaction which 
is more stronger between heavier matter may play a role 
similar to elastic network in producing the spongelike 
large-scale structure. This explanation seems to be con- 
sistent with a standard picture of the universe evolution 
(a gravitational-instability model) that such a heteroge- 
neous structure develops by gravitational amplification 
of density fluctuations. 

This universal appearance of sponge structures in 
phase separation of these systems originates from volume 
phase transition, or more strictly elastic phase separation 
of a dynamically asymmetric mixture that is composed of 
a network-forming component and a fluid ( such as a liq- 
uid and a gas). The elastic network can be a real one 
as in gels (permanent network) and polymer solutions 
(transient network) or a virtual one due to attractive 
interactions. In the former, the real structure having 
large internal degrees of freedom can store the elastic en- 
ergy for bulk compression, while in the latter the virtual 
network due to attractive interactions can also store the 
elastic energy. In this sense, we can conclude that the ex- 
istence of both the component having bulk (relaxation) 
modulus and the fluid component is a prerequisite for 
the formation of spongelike structure due to the volume 
shrinking of one phase. 

For example, such phenomena are never observed in 
solid mixtures, except for the case that the mobility is 
strongly dependent upon the composition Phase 
separation of elastic solid mixtures {e.g. metal alloys) 
does not accompany a drastic volume change of each 
phase if there is no strong composition dependence of 
mobility. This difference causes a marked contrast be- 
tween elastic phase separation in solid mixtures and vis- 
coelastic phase separation; in the former, a softer phase 
always forms a continuous phase to minimize the total 
elastic energy |Q, in contrast to the latter. 

B. Physical origin of volume shrinking 

We briefly discuss the physical meaning of the above 
criteria of the formation of sponge structure, or relative 
volume shrinking of a more viscoelastic phase. This is 
related to the fact that in a two-fluid model V • Vk needs 
not to be zero even under the incompressible condition 



V • w = for the average velocity. There are three impor- 
tant factors in this problem: (i) whether the component 
k is compressible in a mixture, or not, (ii) whether V • vi 
is large enough or not, and (iii) whether the change in 
\7 • ilk is properly coupled with the stress, or not. 

The condition (i) is usually satisfied since we can 
change the spatial configuration of one component ar- 
bitrarily in a two-component mixture, in general. The 
condition (ii) is satisfied only for the system containing 
a fluid as its component. Finally, the condition (iii) is 
satisfied only when there exist attractive interactions be- 
tween the components. In a simple fluid mixture, for ex- 
ample, V • Ufc is not coupled with the elastic stress even if 
there is a difference in viscosity between the two compo- 
nents. In relation to this problem, we consider the case of 
a mixture whose components have different glass transi- 
tion temperatures as an example. If only one component 
becomes viscoelastic, the deformation of this component 

(k) ^ ^ 

that is described by kIj' and V • Vk causes the elastic 
stress. Thus, the more viscoelastic phase becomes the 
matrix phase and forms the spongelike structure as ex- 
perimentally observed Q| . To prevent usual spinodal de- 
composition from taking place, the bulk modulus should 
be sufficiently large: The initial growth of the concentra- 
tion fluctuation has to be suppressed mainly by V ■ Vk. 
For example, this is realized by the formation of transient 
interaction network. 

We believe that the bulk mechanical relaxation mod- 
ulus K{t) plays an essential role in the fluctuation sup- 
pression and the volume shrinking, while the bulk os- 
motic modulus Kos does not play a primary role. The 
interactions can have any origins including entropic and 
energetic ones. The long-range nature of the elastic in- 
teraction in a more viscoelastic phase is directly related 
to how efficiently the fluctuation is suppressed and the 
volume of the relevant phase can be changed and, thus, 
to the ability of the formation of spongelike structure. 

In relation to the above, it should be mentioned that 
the strong composition dependence of mobility has a sim- 
ilar effect since it slows down the diffusion process 
in the slow-component-rich region selectively. Thus, the 
similar behavior of phase inversion has been observed 
even in the framework of a solid model (model B) [ p9[ . 
As discussed in Sec. V. B, the physical factor responsi- 
ble for the asymmetric stress division in a fluid mixture 
is dynamic asymmetry between the two components of a 
mixture. This strongly indicates that the most essential 
physical origin of volume shrinking behavior and the re- 
sulting phase inversion is the coexistence of "asymmetry 
in mobility between the two components of a mixture" 
and "attractive interactions" , irrespective of whether a 
mixture is solid or fluid. 
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C. Roles of shear relaxation modulus on the 
formation of a networklike structure 



We have already discussed the roles of bulk relaxation 
modulus, namely, the suppression of the homogeneous 
growth of concentration fluctuations. Here we focus our 
attention on the roles of shear relaxation modulus. An 
important fact is that the bulk relaxation modulus is 
closely related to the diffusion while the shear relaxation 
modulus is not: The bulk stress gradient V ■ ctb, is usu- 
ally (at least in the initial stage) in the same orienta- 
tion with the osmotic stress gradient, V • 11, since both 
are related to the diagonal part of the deformation ve- 
locity, V described before. On the other hand, 
the shear stress gradient, V • erg, is usually not in the 
same orientation with V • 11, since it is related to the 
off-diagonal part of the deformation velocity. Thus, we 
think that the shear relaxation modulus plays a domi- 
nant role in the formation of a networklike structure in 
the intermediate stage of viscoelastic phase separation: 
The overlapping of the stress fields having the spherical 
symmetry around a spherical solvent holes induce the de- 
formation of shear type. This initial spherical symmetry 
of the stress field is characteristic of the bulk stress fields 
coupled with V • tv. The shear deformation causes the 
shear stress fields through the shear relaxation modu- 
lus. Thus, the thin part of a more viscoelastic phase can 
support the shear stress and be elongated further. In 
other words, the existence of shear relaxation modulus 
is responsible for the formation of a networklike pattern 
composed of highly elongated thin structures |2|]. We 
believe that without shear relaxation modulus, the net- 
worklike pattern with threadlike structures can never be 
formed. 



D. Difference in elastic effects between solid and 
fluid systems 

Here we discuss the difference in phase-separation mor- 
phology between elastically asymmetric solid mixtures 
psf and dynamically asymmetric fluid mixtures. 

In the diffusion-dominated process, the system ap- 
proaches to the final equilibrium state to reduce the total 
free energy including the elastic energy [see Eq. (|30|)]. As 
a result, the morphology that minimizes the elastic en- 
ergy is selected. This is the case of solid mixtures having 
only elastic asymmetry but no dynamic asymmetry. In 
relation to this, it should be noted that the solid mixture 
having dynamic asymmetry behaves entirely differently 
(see e.g. Ref. |3|]). 

In the flow-dominated process, on the other hand, the 
force balance condition plays an essential role in pattern 
selection [see Eq. (^3|)]. As a result, the morphology it- 
self is determined by the force balance condition. The 



asymmetric stress division leads to the spongelike struc- 
ture where the more viscoelastic phase forms a continu- 
ous networklike structure to support the stress. Further, 
the two-fluid nature makes the volume change of phases 
possible. 

Thus, we can say that dynamic asymmetry is a pre- 
requisite to the phase inversion, irrespective of whether 
material is solid or fluid. 



X. CONCLUSION 

In summary, we obtain a general model of viscoelastic 
phase separation on the basis of a two-fluid model: We 
demonstrate that the bulk relaxation modulus plays an 
important role in viscoelastic phase separation even in 
polymer solutions. Our recent simulations based on the 
viscoelastic model indicates the importance of this bulk 
mode [^. Inclusion of this effect makes a viscoelastic 
model quite general: The viscoelastic model can describe 
phase separation or critical phenomena in any isotropic 
condensed matter without any exception, if there is no 
coupling with additional order parameter. 

Although the viscoelastic model is a quite general 
model of critical phenomena and phase separation, the 
critical behavior of this model may be intrinsically 
nonuniversal in the sense that internal slow modes of ma- 
terial can affect the critical dynamics even near the crit- 
ical point. This problem needs further studies to check 
whether the dynamic universality practically breaks in a 
mixture having strong dynamic asymmetry between its 
components or not [^. 

As a straightforward extension of the stress division 



in polymer mixtures |12|, we also propose a simple re- 
lation describing how the stress is divided by the two 
components, on the basis of the idea that the mechanical 
coupling between a component and the mean-fleld rhe- 
ological environment is only due to the friction between 
them. We also point out that the relation is not useful for 
a mixture whose components have different mechanisms 
of the molecular motion. 

We also show that the characteristic features of vis- 
coelastic phase separation can be well explained by the 
concept of " order-parameter switching" between compo- 
sition and deformation tensor. 

We discuss the universal features of spongelike struc- 
tures observed in various material and demonstrate that 
there is a common physical origin that is explained by the 
framework of our viscoelastic model of phase separation. 
It is concluded that the most essential physical origin 
of volume shrinking behavior and the resulting sponge 
structure is the coexistence of "asymmetry in mobility 
between the two components of a mixture" and "attrac- 
tive interactions" , irrespective of whether a mixture is 
solid or fluid. In relation to this, we would like to point 
out that the phase inversion is also observed in a recent 
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simulation of model B including the strong composition 
dependence of the mobility We believe that this is 
only the way to introduce the dynamic asymmetry into 
a solid model. The relation between our model and their 
model and the underlying physics of their similar behav- 
ior [^,|9j will be discussed in detail elsewhere. 
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APPENDIX A: COMMENTS ON APPLICATIONS 
OF VISCOELASTIC PHASE SEPARATION IN 
MATERIAL SCIENCE 



We briefly discuss the application of a spongelike mor- 
phology observed in viscoelastic phase separation. Al- 
though a spongelike structure appear only transiently in 
viscoelastic phase separation, this structure can be frozen 
by suitable methods: (i) simultaneous evaporation of a 
solvent for a polymer solution during phase separation, 
(ii) a further quench of a system below Tg, and (iii) com- 
bination of other processes such as crosslinking reaction. 
We think that spongelike structures reported in litera- 
ture [|5|-|^j40[|| are induced primarily by the mecha- 
nism described here. In relation to this, we point out 
that some sponge phases have periodic structures (see 
e.g. Ref. while others do not as in our case. This 

can be explained by the way of nucleation of solvent holes: 
only when nucleation is heterogeneously induced with a 
high density in a short period, a periodic sponge struc- 
ture can be formed by the long-range elastic interaction 
between solvent holes (correlated nucleation). 

We also point out that polymerization-induced phase 
separation may lead to the networklike structure of a mi- 
nority phase if there is a certain degree of dynamic asym- 
metry induced by the polymerization of a component. 

In the common sense view of conventional phase sepa- 
ration, a minority phase never forms a continuous phase 
and forms only an isolated phase However, our 

present study indicates the possibility that we can in- 
tentionally form a spongelike continuous structure of the 
minority phase of a more viscoelastic phase for any dy- 
namically asymmetric mixture. 
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